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Abstract 



The reheating stage in post-inflationary cosmologies is reanalyzed. New 
techniques from non-equilibrium quantum field theory allow a consistent 
derivation of the equation of motion including the non-linearity of the dy- 
namics. These offer a rationale for the elementary theory of reheating based 
on single particle decay which is seen to be valid only in the linear regime of 
coherent oscillations of the scalar field. 

A new non-perturbative mechanism of induced amplification of quantum 
fluctuations is introduced and studied in detail, both analytically and numer- 
ically. This is a non-linear mechanism that is typically a far more efficient 
way of transfering energy out of the zero mode and into production of lighter 
particles than single particle decay. Thermalization is discussed and we esti- 
mate the reheating temperature to be of the order of the infiaton mass, thus 
providing a potential solution to the Polonyi and moduli problems. 
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I. INTRODUCTION 



In typical inflationary scenarios of the early universe, after the many e-folds of inflation 
necessary to solve the horizon and homogeneity problems, the matter and radiation energy 
density had been red-shifted to almost zero. However, once inflation has ended, this scenario 
has to merge with the standard big bang, radiation dominated cosmology. The succesful 
intertwining of inflation with the standard big bang cosmology thus necessitates some source 
of energy and entropy to rethermalize, or reheat the universe 

Models that reheat to temperatures above the GUT scale would have to cope with the 
monopole and domain wall and inhomogeneity problems associated with the breakdown of 
their symmetries. However there is a fundamental scale for a reheating temperature: nucle- 
osynthesis, this requires reheating temperatures in excess of 10-100 Mev. If baryogenesis can 
occur at the electroweak scale then it would be imperative to have reheating temperatures 
larger than say a few hundred Gev (a baryon over antibaryon excess may be wiped out if 
this occurs during the fast inflationary era). Therefore, nucleosynthesis and particle physics 
restrict the reheating temperature but leave a wide window open for the acceptable values 
of the reheating temperatures: certainly larger than about 100 Mev, and desirably, larger 
than say about 1 Tev but not close to GUT scale. 

The standard picture [0-0] of old and new inflation invokes a scalar field, the inflaton, that 
produces a phase transition (either second or first order) and that at the end of inflation the 
expectation value of its zero momentum mode oscillates about the minimum of its (effective) 
potential. In chaotic inflation |Q, the inflaton energy density still drives an inflationary 
phase, but there need not be any phase transition for this to occur. This inflaton field is 
constrained to have a mass a few orders of magnitude smaller than the Planck mass, 
but typically much larger than the masses of the particles involved in the particle physics 
models. In the standard view of inflationary models, the expectation value of the inflaton 
field oscillates around the minimum of its potential, and its couplings to the lighter particles 
allow the inflaton field to decay. This decay process is then supposed to induce a damping 
term in the evolution equation for the inflaton expectation value of the form r0 with T being 
the total decay rate of the inflaton field The standard estimate for the reheating 

temperature based on single-particle decay 0] is then obtained by comparing the total decay 



rate of the particle (F) to that of the expansion, obtaining Tr ~ 0.1 ^FM^; 0. 

Recent investigations of the non-linear quantum dynamics of scalar fields reveal a variety 
of new and striking phenomena [§-0. The main relevant implication for the reheating 
problem being that the particle production induced by the time evolution of the inflaton is 
significatively different from linear estimates. 

The non-linear (quantum) effects lead to an extremely effective dissipational dynam- 
ics and particle production even in the simplest self- interacting scalar field theory [1TT1,13 



m 



which single particle decay is kinematically forbidden. In the case in which only the classical 
evolution of the expectation value is taken into account in the evolution of the quantum fluc- 
tuations, this mechanism corresponds to parametric amplification or parametric resonance, 
because the classical time evolution is periodic in time [jT5|. However when back-reaction 
effects are taken into account, particle production damps out the evolution of the scalar 
field, and this damped evolution is incorporated in the equations for the quantum fluctua- 
tions. In this case the time evolution of the scalar field including the dissipative effects is not 
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periodic; this implies that the situation is far different from parametric amphfication. We 
refer to this situation, which includes the back-reaction effects as induced amplification, 
to distinguish it from parametric resonant amplification [|T5|. 

Back-reaction effects [Tl|Jl3[] drastically change the picture of particle production. In the 
case of parametric amplification, particle production never shuts off and the total number of 
particles created is infinite in Minkowski. It is only when the back-reaction is incorporated 
that the damping effects on the expectation value feed back to the particle production 
mechanism, eventually shutting it off. This has been studied in detail analytically and 
numerically in |]T3[ for a self-interacting scalar field. 



Induced amplification is most effective when the amplitude of the expectation value of the 
scalar field (from the equilibrium value) is large and is, therefore, an intrinsically non-linear 
process, which we refer to as non-linear relaxation. 

Particle production via parametric amplification of quantum fluctuations has been stud- 
ied recently in the semiclassical approximation, but without back-reaction effects, in con- 



nection with fermion production by (pseudo) Nambu-Goldstone bosons |]T6 



Induced amplification and non-linear relaxation will be the primary mechanism for dis- 
sipation via particle production in the case of large amplitude of the scalar field. Such is 
the case in chaotic infation scenarios |^, as well as in the out of equilibrium regime during 
phase transitions for fields with very fiat potentials near the maximum, which is typically 
required for a long inflationary phase and also occurs for the moduli fields in string theory 
ir?! . In these inflationary models with light moduli fields, with masses typically of the order 
of a Tev, and couplings of order (Tev/Mpi) the standard reheating theory predicts reheating 
temperatures of the order of O.OlMev, well below nucleosynthesis temperatures, this is the 



Polonyi problem |17H20[ . 

The non-linear mechanism of induced amplification of fluctuations and copious particle 
production in the large amplitude regime of the scalar field will necessarily modify the 
standard picture of reheating. The standard reheating scenario was primarily based on 
the premise that particle production only ocurred when the inflaton oscillates with small 
amplitude at the bottom of the (effective) potential. Understanding the time scales and 
non-linear processes of reheating and eventually thermalization acquires further importance 
with the (speculative) possibility that asymptotic oscillations of the inflaton field around 
the minima of the (effective) potential may still be present in today's universe in the form 
of dark matter 

Thermalization is a process that is fundamentally different from particle production. 
Typically, as will be seen below, the particles produced by the process of induced amplifi- 
cation will be in a non-thermal distribution [|T^]. Thermalization is a collisional process, in 
which the produced particles will exchange energy (and momentum) and eventually achieve 
a thermal distribution. In principle the time scales for dissipation via particle production 
may be widely different from the time scales of collisional thermalization. For very weakly 
coupled theories, if the distribution of produced particles via induced amplification is very 
far from thermal, many collisions will be necessary to thermalize the system and the time 
scale for thermalization may very well be large compared to the time scale of dissipation via 
particle production. This will have implications in cosmology discussed in the conclusions. 

The process of thermalization, reheating and relaxation of perturbations, and the ensuing 
production of entropy is also very relevant in heavy ion collisions and in phase transitions 
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in particle physics, at the electroweak scale within the context of baryogenesis [E^ , as well 



as for the quark-gluon plasma, deconfining and chiral phase transitions P2| , ^ . 

In the linear relaxation case the important scale is determined by the decay rate F, usually 
referred to as the damping rate. The literature usually identifies F ^ ImS(ti;, k)/uj with 



the "thermalization rate". Here lmJ]{uj,k) is the imaginary part of the self-energy. In this 
article we offer a real time critique of this relation by studying the real time evolution and 
relaxation of linearized perturbations, and point out the following observations: 

1. The damping rate is identical to the imaginary part of the self energy only when the 
inflaton (the particle interpolated by the order parameter) is a resonance with F being 
its width. In the complex frequency plane this corresponds to a pole in the second 
(unphysical) Riemann sheet. Even in this case, this damping rate corresponds to the 
relaxation of the expectation value of the scalar field and is exponential for some time 
regime. However, eventually relaxation continues with a power law tail. When the 
imaginary part is zero on-shell and the one particle pole is below the multiparticle 
thresholds, relaxation is given by a power law and no damping rate can be associated 
with the imaginary part of the self energy. Moreover, the relation F ^ Im E(ci;, fc)/u; 
only holds in the linear approximation around equilibrium. 

2. Even when relaxation of the expectation value of the scalar field is exponential, this 
damping rate determines the approach to equilibrium of the expectation value of the 
scalar field but it cannot be immediately inferred that the same time scale describes 
thermalization, i.e. the approach to a Bose Einstein distribution of an initial off- 
equilibrium distribution. Thermalization is a rather different process and has to be 
described with a Boltzmann equation with a collision term. Thus, rather than inter- 
preting this time scale as a thermalization scale, we interpret it as a relaxation scale 
for the expectation value. 

3. We consider the non-linear quantum field evolution (non-linear relaxation) incorpo- 
rating the order parameter into the effective mass in a self-consistent way. [The self- 
consistency requirement makes the evolution equations non-linear]. We study the 
inflaton evolution in real time coupled both to light scalars and fermions, providing 
exhaustive numerical results. We find that the non-linear relaxation time scales can 
be much shorter than those predicted by linear relaxation. In addition, the particles 
are produced with a momentum distribution that is very far from thermal and skewed 
towards low momentum. In the case of fermions we provide numerical evidence for 
the phenomenon of Pauli blocking that hinders dissipation and production of fermion- 
antifermion pairs. This phenomenon is very similar to that found by Kluger and 
collaborators [EH] in their studies of fermionic back-reaction in the presence of strong 



electric fields. 

4. We find that the reheating temperature is in between the inflaton mass and the lighter 
scalar mass. The precise relation is argued to depend on the product of the coupling 
and the initial amplitude of the inflaton zero mode and also logarithmically on the 
coupling. The distribution of produced particles is very far from thermal and skewed 
towards momenta of the order of or smaller than the inflaton mass. We argue that 
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this implies fairly long thermalization times and a wide separation of the relevant time 
scales. 



Both time scales will have to be understood in detail to provide a quantitative and reliable 
estimate of the reheating temperature. 

Although we are ultimately interested in describing reheating and thermalization in 
Friedmann-Robertson- Walker cosmologies, in this article we will work in Minkowski space 
assuming that the relaxation time scales are much shorter than the expansion time scale. 
Furthermore, we want to isolate dissipative effects arising from particle production from 
those resulting from the red-shift in an expanding cosmology. 

A more detailed account of the results presented here can be found in ref. p7 . 



In section II, we introduce the general model that we propose to study and summarize 
the necessary ingredients of non-equilibrium field theory to provide the reader with the 
essentials needed to reproduce our calculations, as these do not seem to be part of the 
standard techniques. In section III, we study linear relaxation in real time and elucidate 
the role of the imaginary parts of self energies, their interpretation in real time and criticize 
the identification with "thermalization rates" . Here we provide a derivation of the reheating 
temperature based on single particle decay, valid in the linear regime of coherent oscillations 
of the scalar (infiaton) field. We also discuss the basics of the Polonyi problem. 

In section IV we study the process of non-linear relaxation, and particle production 
both for scalars and fermions, providing exhaustive numerical results. Here we also discuss 
the issue of thermalization, and the important time scales associated with this collisional 
process. Under the assumption of thermalization rates larger than the expansion rate we 
provide an estimate of the reheating temperature based on the total number of produced 
particles. 

Finally we conclude in section V with a discussion of the implications of our results, in 
particular we emphasize the possibility that the new mechanism for reheating may provide 
a solution to the Polonyi problem. An Appendix is devoted to a pedagogical exercise in 
non-equilibrium field theory. 



II. THE MODEL AND THE METHODS 

We consider the simplest model [QJ^ where the infiaton field $ couples to a scalar a and 
to a fermion field ip. That is, 

C = -U{^' + ml+ga')<t>-^<l>'-^cx{^' + ml)a-^a^ + ^P{^^-m^-y<^)^ . 

(2.1) 

The case ma, <^ m$ will be of particular relevance, since in this case there are open 
decay channels for the infiaton. 

We will investigate the scalar and fermionic couplings independently. Although the 
situation for reheating corresponds to (almost) zero temperature, we will study the case of 
linear relaxation at finite temperature. The reason for this is that finite temperature effects 
reflect the Bose enhancement and Pauli blocking factors that appear whenever there are 
excitations in the medium. These contributions from the medium will allow us to identify 
similar physical features in the case of non-linear relaxation. 
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A. Linear Relaxation: Amplitude and Perturbative expansion 



To explore the behavior of the inflaton within the linear regime, we first use the tadpole 
method to obtain the equation of motion (see the Appendix for details on how this 
method is actually applied and reference ||28| for an alternative implementation). The next 
step is to linearize this equation in the inflaton zero mode amplitude and use this to study 
the relaxational dynamics of the inflaton. We should note that this amplitude expansion is a 
priori different from the standard perturbative expansion in the relevant coupling constant. 
Later in this work, we will compare the results obtained here with results obtained through 
a self-consistent, non-perturbative resummation both in the coupling constants and the field 
amplitude. 

We arrange the initial conditions to be such that the fields a, ip start to interact at a 
time that we choose t = 0. This can be accomplished by making the coupling "constants" 
time dependent, i.e. zero for t < and different from zero for t > 0. To perform the 
calculations, we will need the non-equilibrium Green's functions and Feynman rules. Since 
the non-equilibrium generating functionals involve a forward and backward time contour 



i8| j25|j29[| , the number of vertices is doubled. Those in which all the fields are on the forward 
branch (fields labeled by (+)) are the usual interaction vertices, while those in which the 
fields are on the backward branch (fields labeled by (— )) have the opposite sign. The 
combinatoric factors are the same as in usual field theory. The spatial Fourier transform of 
the necessary finite (initial) temperature propagators are: 

• Bosonic Propagators 

G++(t, t') = G>{t, t')e{t - 1') + G<{t, t')e{t' - 1) , 
G,-(t, t') = GUt, t')e{t' -t) + G<{t, t')e{t - 1') , 

G+-(t,t') = -G',<(t,t'), 

G,+M = -G>M, (2.2) 
i 



{[l+n,(a;fc)]e-^-'=(*-*') + n,(a;fc)e*-'=(*-*')} , 



2ujk 

G^{t,t') = i fd^x e-'^-"^ {<l>{6,t')<!>{x,t)) 



2ij 



{[1 +n,(.;,)]e^-'=(*-*') +n,(^,)e-^-'=(*-*')} 



where m is the mass of the boson. 

Fermionic Propagators (Zero chemical potential) 

si^M = s^it,t')eit - 1') + s5(t,t')e(t' - 1) , 
5g-(t,t') = 5|(t,t')e(t' -t) + s^{t,t')eit - 1') , 

St-{t,t') = -SSit,t'), 
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S^Ht,t') = -S^{t,t'), 

S^{t,t') = -i [ d^x e-'^-"^ (V'(x,t)V'(0,t')) 



(2.3) 



2uJk 

Sf{t,t') = t f£x e-'^-^ {^{Q,t')^{x,t)) 



+ m^)(l -nf{ujk)) + e 



7o(/; - m^)-fonf{uk) 



nf{ujk) 



el3^k + 1 ■ 



In the linear amplitude approximation, corresponding to linear relaxation, we find in all 
cases the following form of the equation of motion for the expectation value (see Appendix) 



Kjj{t - t') 6j;{t') dt' = 



(2.4) 



K,it-t') = KAt-t'Mt-t') ^ 



where we have imposed as boundary conditions that the infiaton and the other fields are 
coupled at time t = but uncoupled for previous times, and introduced 



d X e 



ip-x 



(x, t) , = p ^ + m| + 5m{T) , 



(2.5) 



where 6m{T) is the time independent (but temperature dependent) contribution from tad- 
pole diagrams. These contributions renormalize the mass and introduce a temperature 
dependent effective mass and will be specified later in each particular case. The quantity 
^r,p{t — t')Q{t — t') is the retarded self-energy. It will be computed to dominant order in the 
couplings for both fermions and bosons. Although the self-energy at finite temperature has 



been computed before in the literature pO| , |3T| , we differ from previous treatments in that 
we perform our calculations directly in real time, this allows us to study real time relaxation 
as an initial condition problem. 

Eq. ( p.4|) can be solved by Laplace transform. Define 



~st 



6p{t) dt . 



Then, eq. (|2.4| ) becomes 



s Vp(s) - S(5p-(0) - 4-(0) + ^^p- V5p(s) + V5p-(s) Sp>(s) = 



(2.6) 



(2.7) 



with Sp'(s) the Laplace transform of Er^p-(t). 

For computational purposes, it can be shown that at zero temperature each graph of Sp'(s) 
is exactly equal to the corresponding graph of the zero-temperature equilibrium Euclidean 
quantum field theory, s being the time component of the Euclidean four momentum. 

In general Sp-(s) can be written as a dispersion integral in terms of the spectral density 
p{Po,P,T) 
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= ~J ^T^^ "^Po . (2.8) 

The imaginary part of the self-energy is found to be 
ImSp'(s = iuj ± 0^) = ±S/p-(co') , 

j:ip{uj) = n sign {iu)[p{\iu\,p,T) - p{-\u\,p,T)] . (2.9) 

The presence of sign {u) in the above expression characterizes the retarded self-energy. 
Let us choose 5p(0) = 6i , S^{0) = for simplicity. We get from eq. (|2.7] ) 



The Laplace transform can be inverted through the formula 

ioo+e , ^ ds 



5,it) = I e^* if.is) — , (2.11) 



loo+e 



where e is a positive real constant (Bromwich countour). Thus we need to understand the 
singularities of [s^ + Vt^+ Sp-(s)]~^ We now consider the following cases: 

• The infiaton potential admits spontaneous symmetry breaking (SSB), and is only 
coupled to lighter scalar fields. 

• The infiaton is coupled to fermions only. 

We will also consider the subcases where the initial temperature is taken to be zero, as 
would be appropriate in the case of evolution in the post inflationary universe, as well as the 
situation where the initial temperature is non-zero, which would be relevant to the situation 
of a scalar field starting in an initial (non-equilibrium) but thermal state and evolving out 
of it. 



B. SSB with Coupling to Lighter Scalars Only 



If m| = — /i^ < 0, then the new minimum is at $o = y^/^^/"^*' ^"^^ write $^(a;, t) = 
$0 + 0(a^7 i) + X^(a^, i)- The masses are now shifted to 

M^ = 2fj,\ Ml = ml + g'^^l . (2.12) 



The contribution from the quartic infiaton self-coupling has been studied previously |T3 
thus it will not be repeated here. 

The tadpole correction to the mass in eq. ( p. 51) is given by 

d^k , f d^k l + 2nb{ujk) 



SM{T) = -zg J --y^Gt;At,t) = 9 J 



(27r)3 ' ' (27r)3 2uJk 



LOu = \^k^ + Ml. 
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This mass renormalization does not influence the dynamics. The retarded self energy is 
found to be at order g^, 



Kj^{t - t') = 2z^2$2 



(2vr) 



- '')GU,Jt - t') - Glit - t')G5^.^(t - t')\ . (2.13) 

With the non-equilibrium Green's functions defined above, we find 

d^k 1 



(27r)3 2uj^,uj^,^. 



—2 nh{uJk) sin[(co' 



k ^k+p 



{[1 + 2 UbiuJk)] sin[K + ^A?+p)(i - t')] 
){t-t')]}. (2.14) 



The Laplace transform can be written as a dispersion integral in terms of the bosonic spectral 
density PbiPoyP^T) (see eq. (^.8])). We have to one loop level: 







d^k 



d^k' 



{2'Kf2uJkJ (27r)32cufc' 



{2TTf5^{p- k - k') 



X [5{po - t^fc - uJk'){l + 2nb{uJk)) - 5{jpo -ujk + ujki)2nh{ujk)] ■ 

The imaginary part of the self-energy is given by eq. ( p.9| ). 

We will analyze only the case of a spatially constant order parameter corresponding to 
p = because we will later compare with the case of non-linear relaxation which we only 
study for a homogeneous (translational invariant) expectation value. In this case the spectral 
density can be written as 



PbiPo,0,T) 



l + 2nj^ 



PbiPo, 0,0) . 



(2.15) 



The spectral density pb{po, 0, 0) is a Lorentz scalar and is proportional to the decay rate of 
the boson $ into two a particles. It is a straightforward exercise to find 



Pb{Po, 0,0) 



2^2 




1 - 



AMI 

pI 



e [pt - AMi 



(2.16) 



Clearly Sq(s) has a logarithmic divergence, independent of s and T. We choose to subtract 
this divergence at s = and absorb the subtraction into a further temperature dependent 
renormalization of M. In order not to clutter notation, from now on we will refer to M as 
the fully renormalized mass including the above subtraction, and to E(s) = S(j(s) — Sg(0). 



From equations ( p.8| , p.l5| , |2.16| ) we find that the self-energy has an imaginary part above 
the two cr-particles threshold given by eq. ( p.9|) with 



Moo) 



9'% 



57r 



1 - 



AMI 



e cj^ -4M^ sign (cj 



(2.17) 



This expression is recognized as the imaginary part of the retarded self-energy (determined 
by the sign (a;)) and shows the usual Bose enhancement factor 
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1. Zero Temperature 



In this case from [2.141 we find 



4vr2 



4M2 , , 1 
1 + ArgTh 



1 + ^ 



(2.18) 



In order to compute the inverse Laplace transform through eq. ( |2.11| ) we must first study 
the analytic structure of ip{s) in the s-plane. (p{s) has poles at the zeroes of 



s^ + M^ + Eis) = 



(2.19) 



These correspond to $-one-particle states with the mass including one-loop radiative cor- 
rections. At zeroth order the poles are purely imaginary 



s± = ±iM . 
To find the one-loop correction, we set 

s+ = iM + r 

and similarly for s_. Inserting this in eq. ( p.l9| ) yields to order g"^ 

2iMr + T.{iM) = . 

That is, 

.S(zM) 



r = i- 



2M 



(2.20) 



f2.211 



f2.22) 



(2.23) 



When M < 2M„, S(iM) is real (see eq. (|2.18| )) and eq. ( p.23|) gives a real correction to the 
$ mass 



s± = ±iMo = ±i 



M + 



ArgTh 



1 



4M2 
M2 



- 1 



(2.24) 



The Laplace transform ip[s) also exibits a cut along the imaginary axis starting a.t s = iu 
±2iMo-. For s in the first Riemann sheet (physical sheet) we obtain 



Sphysical(^tU ± 0+) = T.r{uj) ± iT^i^u) 

with Hpi and S/ both real and given by 



(2.25) 



Sr(u;) 



47r2 



1 



4M2 



^ ArgTh Wl 



4M2 



S/(^) 



4M2 
1 ^ > . 



(2.26) 



We can now proceed to compute the inverse Laplace transform ( |2.11| ) by deforming the 
contour. 
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m 



Si COS Mot _^ 2Si 



TC 



2 A/a \0J 



(2.27) 



For large time M^t ^ 1 the integral over the cut is dominated by the endpoint u = 2Mfj 
and goes to zero as 

Scut (t) 



6i M2 g^<!>l cos(2M,t + 



47r2(M2-4M2 + ^)2 



(M.t)i 



f2.28) 



The is completely determined by the behavior of the spectral density at threshold. 

The situation changes drastically for M > 2M„. In such case the ^-particle is unstable 
and thus S(iM) becomes complex and its value depends from which side of the cut we 
approach the imaginary axis. Now the solution for the pole will be complex and r will 
acquire a real part. Due to the discontinuity in S/ across the two-particle cut, eq. (p.23|) 
must be written carefully as 



2M 



S(iM + Re[r]) 



(2.29) 



where the (small) real correction inside the argument will determine on which side of the 
cut the solution resides. Then the real part of r should satisfy 



Rer 



1 



2M 



ImS(iM + 0+ sign (Re [r])) . 



(2.30) 



From eq. ( |2.171 ) we see that ImEphysicai is negative for sign (Re [r]) < and positive for 
sign (Re [r]) > 0. Therefore eq. ( p. 30 ) has no solution in the physical sheet. 



The analytic continuation of S into the second Riemann sheet is such that 



E^^iu ± 0+) = ^r{u) t i^iiuj) 



u > 2M„ 



and we find the solution 

1 



Rer 



Imr 



2M 



2M 



levrM 



AMI 



<0, 



Stt^M 



AMI 



ArgTh \1 



AM^ 



(2.31) 



For M > M„, 



Imr 



Sti^M \^M„ 



(2.32) 



I Re [r]| coincides with the decay rate $ — > 2(T (as it must be) which is the rate per unit time 
to produce a particles. The poles s± move off into the second sheet when M becomes larger 
than 2M„ as expected |3^ . 



Thus when we compute the inverse Laplace transform ( |2.11| ) in the unstable case we are 
left with the integral over the cuts since both poles are in the second Riemann sheet and we 
find 
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m 



2(5,; 



ujUjIuj) cosut duo 



TT J2M. [UJ^ - M2 - Sr(cu)]2 + 



(2.33) 



Since now M is inside the integration region, for weak coupling there is a narrow resonance 
at C(j ~ M. Thus for weak couphng it takes Breit-Wigner form and we find to a very good 
approximation, 



5{t) ^5iA e-^*/2 cos(Mt + a) 



where 

A = l + 
For M > M^, 

A = l + - 



dT.R{M) 



SvrM 



4M2 
1^ ' 



9E/(M) 



(2.34) 



(2.35) 



M 



a 



+ 



M 



(2.36) 



The Breit-Wigner approximation, however, is vahd only for times < F"^ ln(F/MCT); for longer 
times the fall off is with a power law t~^/2 determined by the spectral density at threshold 
as before. 



2. Non-Zero Temperature 



The physical mass gets a finite temperature dependent correction from the tadpole con- 
tribution 6M{T) given by 



6M{T) = g r 
Jo 



dk 



fc2 



1 



/o 27r2 ^^2 + Ml 6^^^='+*^' - 1 
For small /3 (T ^ M^) this correction takes the form 



(2.37) 



5M{T) = g 



T 
12 



4tt 8tt Mrr 



(2.38) 



We will assume that gT^ <^ M'^ so that we are in the perturbative regime, since otherwise 
hard thermal loops must be resummed ||3^ , |35[| , a task beyond the scope of this article. 

The imaginary part of the self-energy is given in eq. (|2.17|) and the real part can be 
obtained from the dispersion integral ( p.8| ) using ( |2.15| , p.l6| ). The real part of the self-energy 
is difficult to compute for arbitrary temperature, and we just quote its large T behavior: 



\ 



1 + 



AMI 



+ C(T°) . 



(2.39) 



It is interesting to compute S(t,/3) in configuration space for large T. We have by inverse 
Laplace transform 



12 



/+ioo+e rjq 
e^*S(.,/3) — . 

Upon contour deformation we find 

E(t,/5) = --^ / — y^^^ sin(2M,xt) . 

- Jl X 



(2.40) 



(2.41) 



This function is obviously not concentrated at t = 0. It can be related to a Ji Bessel 
function as 



S(t,/3) 



9^ 

271 



AM /■* /f \ 

1 t-2M„ dx - - 1 Ji(2 M^x 



(2.42) 



We find for small and for large t 



2n 



^ 4M^ 2 
1 -t + 0{f 

TT 



g^%T sin(2M^t-7r/4) 



27r2 



(M^t)3/2 



l + 0(t 



The behaviour of the kernel S(t, (3) shows that it cannot be approximated by a phenomeno- 



logical term V ^ even for high T. (As shown in ref. |T3[ this is not the case either for T = 0) 



We can now repeat the analysis of the previous section for both cases M < 2M^ and 
M > 2M(j. In the first case, the one particle pole is below the two cr-particles threshold with 
a (small) finite temperature correction since we are restricted to the perturbative regime in 
which gT^ <^ M. The long time behavior will be oscillatory with the frequency correspond- 
ing to the one particle pole plus long-time power law tails similar to the zero temperature 
case. The second case is more interesting, since now the scalar "infiaton" is unstable, and 
the pole moves off into the second Riemann sheet. In the physical sheet there is a resonance 
with a finite temperature width given by 



r(T) = r(o) 



m 



8nM 



'1- 



4M2 



(2.43) 



The Bose enhancement factor increases the rate and therefore enhances dissipation via the 
production of particles. Although this factor arises from the thermal distribution, we expect 
in general that whenever there are bosonic excitations present, the relaxation rate will be 
enhanced as a consequence of bose statistics, independently of whether these excitations are 
thermally distributed. 



C. Unbroken symmetry with coupling to light scalars 

In the unbroken symmetry case M<j, is above all thresholds and hence $ is a stable 
particle. The order parameter is again given by a formula like eq. ( |2.27|) except that now 
the integration starts at u; = M + 2M„. 

The first perturbative contribution to the kernel S(t — t') in the infiaton equation of 
motion (|2.4|) are now the two loop order graphs usually called "setting sun". Since these 
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graphs are quite complicated, we only perform the two loop computation at zero temper- 
ature, for which we can exploit the relationship with usual Euclidean field theory. Even 
with this simplification the computation is hard, and we limit ourselves to the evaluation of 
the imaginary part of the retarded self-energy near the branch point. We have seen below 
eq. ( |2.27| ) how this actually determines the long time behaviour of the order parameter. In 



general, if T,i{uj — > ciJthrcshoid) vanishes as (u; — u^threshoid)", then 5cut(^) decays as t^-*^^" for 
large times. 

We find for the two loops self-energy 

EK. ^ M + 2M„) . ^^^^J^^^l^^ - (M + 2M„)T . (2.44) 

This yields for Scy^tit) a power decay for long times. 
For massless a we find 

2 

yielding a 5cut(^) decaying as for long times. In summary, in the unbroken symmetry 
case 5{t) is given by the one-particle term oscillating with frequency Mq plus the power 
damped cut contribution (5cut(^)- 



D. Inflaton coupled to Fermions 



Fermions can be treated similarly in the broken and unbroken symmetry phase. To treat 
both cases on equal footing we now define M as mass of the inflaton (scalar) field and m as 
the fermion mass in either case. Using the Feynman rules described in the first section, we 
obtain to one-loop level the kernel 



- t') = ty^ 
^r,p{t - t') = iy^ 



(2^ 



Tr 
Tr 



zS^it, t')^SSJt', t) - zSS{t, t')tSl.{t', t) 



We now concentrate on the homogeneous case p = 0. Using the fermionic Green's functions 
given in the first section it is straightforward to find the fermionic spectral density to be 
used in eq.( p.8D 



PfiPo,T) 



o 9 Po 



Am 



[pi - Am' 



From eq.( |2.9|) , the imaginary part of the self energy is then 



Sir 



1 - 



Am 



1 - 2n/ ( - 



- Am' 



The finite temperature factor reflects the Pauli blocking term |^. It is clear that the zero 
temperature part of S(s) diverges quadratically and that two subtractions are needed. The 
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first one is independent of s and contributes a (quadratically divergent) mass renormaliza- 
tion. The second one is logarithmic divergent and consists of an s independent term that 
adds to the mass renormahzation and another proportional to that will be absorbed in 
wave function renormahzation. We choose to subtract at zero temperature and at an arbi- 
trary scale K. The Laplace transform for the zero momentum component of the equation of 



motion (|2.1CI|) becomes 



Si s 



+ i In + MfniT) +y'll {s, T, k) 



(2.46) 



where Mf^ contains the mass renormahzation and 11 is the twice subtracted kernel, which 
depends on the renormahzation scale. Defining 

MR{T,K) = zj M,n{T) , 

(fR{s,K) = Z^^ lf{s) , 

we finally obtain the renormalized Laplace transform 



Si s 



s' + M^^{T,K)+yU^)U{s,T,^) 



(2.47) 



The inverse Laplace transform is obtained as in equation (|2.11|) . The result will be the 
function Sj^fj{t, k) which is not a renormahzation group invariant. It is clear, however, from 
the renormahzation prescriptions described above that ratios of the amplitude at different 
times, such as S^^{t, k) /Sj^^{Q, n) are renormahzation group invariant. Now the analysis 
can proceed as in the previous section. To obtain the inverse Laplace transform we must 
recognize the sing ularities in (lO^ ). If Mr < 2m (m is the fermion mass in the loop), there 
is a one particle pole (with strength different from one because we decided to renormalize 
off-shell) and a two-fermion cut at = — 4m^. At long times the amplitude will oscillate 
with an oscillation period given by the position of the pole, which is perturbatively close to 
Mji. The contribution of the cut falls-off at long times as and is determined by the 
behavior of the spectral density near the two-fermion threshold. 

More interesting is the case in which Mr > 2m. As in the bosonic case, the pole moves 
off the physical sheet into the second sheet. In the physical sheet the spectral density at 
weak coupling features a sharp peak at Mr + 0{y^) with width 



r(T) = r(o) 



1 - 2n 



M. 



m 



Svr 



-M, 



R 



4m^ 



(2.48) 



A Breit-Wigner approximation predicts exponential relaxation but eventually at long times 
(an estimate similar to the bosonic case) a power law relaxation ensues, completely 
determined by the spectral density at threshold. 
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Pauli blocking makes the resonance narrower and the hfetime of the decaying particle 
longer. The interpretation of this phenomenon is simple. In the thermal bath, fermionic 
excited states are filled with the Fermi-Dirac distribution. In order for the scalar field 
to decay, it must create a fermion-antifermion pair. However, at finite temperature, the 
available states are already filled with thermal excitations and because of the Pauli exclusion 
principle, are not available. At infinite temperature (and zero chemical potential), each 
fermion and antifermion state are populated with occupation 1/2 per spin degree of freedom; 
in this limit the decay rate goes to zero, the lifetime to infinity and the bosonic particle simply 
cannot decay because there are no states available to decay into. Even at zero temperature, 
but in a situation in which excited states are occupied, dissipative processes mediated by the 
production of fermion-antifermion pairs will be hindered by the Pauli exclusion principle, 
since states will already be occupied and no longer available in the particle production 
process. In highly excited states, we expect damping via production of fermion pairs to be 
strongly suppressed by Pauli blocking. This phenomenon has been seen numerically in the 
case of fermion pair production in strong electric fields p6| and will be seen numerically in 
the non-linear relaxation case later. 



E. Thermalization or Relaxation? 

From the analysis presented above, the following conclusions for the linear regime become 
very clear: 

1. The imaginary part of the self-energy only determines a relaxation rate in the case of 
a resonance, that is when there is an imaginary part on-shell for the external particle 
and the (quasi) particle pole moves off the physical sheet into the second (unphysical) 
sheet. In this case a Breit-Wigner approximation describes the exponential relaxation 
for a long time, but eventually the amplitude falls-off with a power law in time, with 
the power determined by the behavior of the spectral density at threshold. When the 
on-shell pole is below the two-particle threshold, the imaginary part does not translate 
to a damping rate. Relaxation is described by a power law with an asymptotic behavior 
completely determined by the position and residue of the pole. 

2. Even in the case of a resonance and exponential damping, the "damping rate" P 
describes exponential relaxation for the expectation value of the scalar field. The issue 
of thermalization is completely different. Thermalization corresponds to the time 
evolution of the (quasi) particle distribution function towards a thermal distribution 
which, in principle, has nothing to do with the relaxation of the expectation value of 
the field. 

An alternative way to look at thermalization is as a process of momentum and energy 
transfer. Thus, the thermalization rate should be identified with the rate of energy and 
momentum transfer which is not necessarily related to the relaxation rate of the expectation 
value of the field. In order to understand thermalization, a collisional (quantum) Boltzmann 
equation must be set up. Although in the Born approximation the collision term includes 
the scattering cross section that is related to the decay rate, the solution to the Boltzmann 
equation implies a resummation that is quite different from the resummation of the Dyson 
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series for the propagator. If the initial distribution is very far from thermal and many 
collisions are necessary for thermalization, the relaxation and thermalization time scales 
may be widely different and in principle unrelated. Thus we insist that the "damping rate" 
obtained from the imaginary part of the self-energy on shell must be interpreted as the 
relaxation rate for the expectation value of the scalar field and in principle not with the 
thermalization rate of the particle distribution function. Moreover, this relaxation rate 
holds in the linear regime (small field amplitude). 



F. The Reheating Temperature: 

We can use the results obtained above to provide an estimate of the reheating tempera- 
ture based on single particle decay and valid only in the linear regime [0,0. Furthermore, 
since the calculations of the relaxation and decay rate were performed in flat Minkowski space 
these will only be useful for an estimate of the reheating temperature whenever T ^ H with 
H being the expansion rate of the universe. Since the calculation of the rate even in flat space 
relies on a "sharp resonance approximation" (Breit Wigner) this implies that the lifetime 
of the particle is much longer than its oscillation period. This translates into the constraint 
rrifp ^ r ^ H. The infiaton oscillates coherently many times around the minimum of the 
potential before its oscillation amplitude diminishes appreciably both as a result of decay 
into lighter particles and because of the red-shift of the energy. 

The Friedman equation for the rate of expansion is 

H^-it) = ^p(t) (2.49) 

with pit) the matter energy density. 

After the reheating stage is completed it is expected that the light produced particles 
will be in the form of radiation (ultrarelativistic), thus p{t) ~ T^eh- During the period in 
which the relaxation rate (or rate of particle production) is F < H{t) the produced particles 
will be red-shifted in energy and will not contribute to a radiation dominated universe. Only 
when the expansion rate becomes smaller than the rate of particle production is when the 
produced particles will remain as a radiation dominated component. Thus the important 
time scale for which the produced particles will remain as a radiation dominated component 
in the universe is that for which H{tf) = F. At this time 

Pitf) = -i^H^itj)^T^\ (2.50) 
on 

leading to the estimate of the reheating temperature within the linear regime: 



- reh 



MpiT (2.51) 



G. The Polonyi problem: 

One area of current interest in which this work may have some implications is that of the 
so-called "post-modern Polonyi problem" concerning fiat directions for some of the moduli 
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fields in string theories [1T7|-[20|]. These are fields with pertubatively fiat directions whose 
degeneracy is lifted by non-perturbative effects. Their masses ^ 10^ — 10^ Gev with 
vacuum expectation values of order 10^^ — 10^^ Gev and couplings to normal matter are 
of order g ^ m^/Mpi ^ 10"^'^. These properties allow the energy density in these fields to 
dominate that of the radiation in the universe until times well after nucleosynthesis. These 
extremely weak couplings determine that the decay rate of these fields is of order 



m5 



r - T7f (2-52) 



A calculation of the reheating temperature based on the "old" scenario of one particle 
decay gives Treh ~ 0.01 — 0.1 Mev , well below the nucleosynthesis scale. 



III. NON-LINEAR RELAXATION 



In this section we study the equation of motion for a homogeneous order parameter 



beyond the linear approximation. This is achieved as follows [|T3[. We write the scalar 
infiaton field as 



$±(f,t) = 0(t)+X^(x,t) 



(3.1) 



with (j){t) the expectation value in the non-equilibrium density matrix and < x^i^^t) >= 0, 
we consider that < >= (this is a consistent assumption). The non-equilibrium path 
integral requires the Lagrangian density 



a 



£(0;X+;a+) = -[(9,x+)'-M2(t) ix^f 



+ 2 



9 



A 

+ 6 



o 



{d,cr+y - m\t) (a+)2j -g<^x^ {o^f - \ ix^Y {^^? + V^+(^/^ - m^{t) - yx^)^' 
A 



M^{t) =ml + -(j)\t) , m\t) = ml + g (j)\t) , m^{t) = + y (f){t) . 

The difference with the linear relaxation case is that we now incorporate 0(t) in the definition 
of the time dependent masses. The equations of motion are obtained as in the previous 
sections, by treating the linear, cubic and quartic terms as perturbations. The necessary 
Green's functions are constructed from the homogeneous solutions of the quadratic forms. 
We again treat the cases where the infiaton is coupled to scalars or fermions separately. 



A. Infiaton Coupled to Scalars Only 



The Green's functions for the scalars are obtained from the mode equations that solve 



dt 



^ + e + M\t) 



Uk{t) = 



(3.2) 



UkiO) = 1 , f/fc(0) = -tWk = -i^k^ + M2(0) 
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Vk{0) 



+ ]? + m^{t) 
1, M^) 



Vk{t) = , 

= —iwk = —i\JkP' + m^(0) 



(3.3) 



The initial conditions on the mode functions (|3.2| , p.3| ) correspond to positive frequency 
solutions at the initial time. In terms of these mode functions, the Green's functions are 

0i 



GUt^t') 



2Wk 
G^^ki^',t) 



[{l+n,{Wk)) Uk{tM{t')+nh{Wk) U;{t)Uk{t')] 



2wk 
GUt\t) 



[(1 + n,{wk)) Vk{t)V:it') + n,{wk) V;*(t)\4(t')] 



and the rest of the Green's functions are given by the relations in equations (|2.2| ). The 
Green's functions above correspond to the situation in which the initial density matrix is in 
equilibrium for the positive and negative frequency modes at time t = 0. This is determined 
by the initial conditions on the mode functions above at an initial temperature T = 1/13. 
However, for the rest of the analysis we will take T = 0. We see that 



(27r)3 2Wk 
d^k |V-fc(t)p 
(27r)3 2wk 



coth 
coth 



2 



Finally the equation of motion to one-loop order for the expectation value is 



A 



A 



(j){t) + mi (j){t) + - <P\t) + - <p{t) < x'ix, t) > +g (j){t) < a\x, t) >= 
D 2 



(3.4) 



If one wants to solve the equation of motion (|3.4| ) to order (h) one would expand in a 
power series in h and only keep the zeroth order term in the mode equations ( |3.2| , ^.3| ). As 
was observed previously [|1^, such an expansion will result in secular terms and becomes 
unreliable at long times. A resummation is necessary to capture the long time behavior 
consistently. We will perform a non-perturbative resummation of the one-loop terms by 
solving the set of equations (^]2| - with the full value of in the mode equations. 
This then becomes a set of coupled non-linear integro-differential equations that provide 
a non-perturbative resummation of select one-loop terms as can be seen by looking at a 
perturbative expansion of the solution to these equations. 

The integration of these coupled set of one-loop equations provide the lowest non- 
perturbativeself-consistent resummation that is energy conserving (see below for energy con- 
servation). This one loop approximation becomes exact in the large (N,M) limit with N 
scalar a fields and M fermionic fields when the a field self-interaction is neglected. 

We want to emphasize this point. By incorporating the full value of (p in the evolution 
equation for the mode functions we are incorporating back-reaction effects self-consistently. 
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If only the classical evolution of is used in the mode equations, we would have parametric 
resonant amplification since the effective frequencies for the mode functions are periodic 
due to the fact that the classical solution is periodic with constant amplitude. This leads 
to particle production that never shuts-off. However particle production leads to damping 
in the evolution of (j). Introducing this damped evolution in the mode equations leads to a 
behavior rather different from parametric resonance: as the evolution of (j) is damped, the 
amplitude becomes smaller and particle production should diminish and eventually stop. 
This was found to be the situation in the case of the scalar field with self-interaction . 



Clearly the integration of this set of equations will have to be done numerically. Because 
the self-interacting scalar case was already studied before |T^, we will now concentrate on 
the interaction with the scalar fields in which thresholds are present. But before we do this 
we must understand the renormalization aspects of this system of equations. 

The large-fc behavior of the mode functions can be understood from a WKB analysis, 
the details of which had already been discussed in [§|,|T3. We find the following divergence 
structure doing this 

- 'ml, + g,(l)\t))\n- + F2{t,K) , 



where A is an upper momentum cut-off and k and arbitrary renormalization scale. The 
subscript h refers to bare quantities and the quantities -^1,2(^5 1^) are finite in the limit A 00. 
We can now read the mass and coupling constant renormalizations 

W^b 2 ^ 2 \ i A 




8^2 12 "'^ '^'V K 



We introduce a further, finite renormalization, by subtracting the functions -Fi,2(^ = 0, k), 
absorbing this subtraction into a (finite) renormalization of the mass. After these renor- 
malizations, we finally arrive at the renormalized set of evolution equations in terms of 
renormalized quantities. We drop the subscript R for renormalized quantities to avoid clut- 
tering the notation, but with the understanding that all quantities are renormalized at the 
scale k: 

2 , , , A 3 A /-A d^k |f/fc(t)P - 1 ^ /""^ d^k \Vk{t)\'^ - 1 



(j)(t) + mi(f)(t) + -(f)-^(t) + - (f)(t) - — - ' ^ " h q Mt) , , 

rv;^ *v^v ; ^ g V ; ^ 2 ^7 (27r)3 IW^ ^ ^ > } (27r)3 2w 

4>{fy\4?{t)-4?{Si)W- = ^ . (3.5) 

To this order, we can replace masses and couplings by their renormalized values in the 
equations for the mode functions (|3.2| , p.3|) . The set of equations ( p.2| , p.3| , |3l5|) provide a self- 




consistent, non-perturbative set of integro-differential equations with back-reaction. This last 
point is extremely important, a periodic evolution for the scalar field in the mode equations 
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( p.2| , |3.3| ) would lead to parametric amplification |Tl| and infinite particle production . 
However the back-reaction of the quantum fluctuations onto the evolution of the zero mode 
of the scalar field is determined completely by these mode equations. When the amplitude 
of the quantum fluctuations grow, this mechanism takes energy of the zero mode, whose 
evolution now becomes damped in amplitude. This damping , a result of the back reaction, 
in turn makes the growth of the quantum fluctuations to diminish. Eventually the particle 
production mechanism must shut-off as a consequence of this back reaction. This is dramat- 
ically different from the catastrophic particle production found by Yoshimura precisely 
because in that analysis back reaction was not taken into account. 

We will now chose the renormalization scale to be k = m$ so as to have only one 
scale in the problem, which makes the numerical evaluation easier. Renormalized non- 



equilibrium equations had also been obtained by Cooper et. al. [^] within the context of 
dynamical evolution during the chiral phase transition. The renormalized equation of motion 
( ^.5| ) may be written without reference to the cutoff A which in the end must be taken to 
infinity, however, numerically the fc-integrals must be calculated with an upper momentum 
cutoff. One must ensure that this numerical cutoff, to be identified with A in the evolution 
equation, be much larger than the masses and amplitudes of the field for the integrals to 
reach their asymptotics and the cutoff dependence to dissapear. We have been careful to 
make exhaustive checks that the final results were insensitive (to the working accuracy) to 
the cutoff which was typically chosen to be A ^ 100 |m$|. This implies that the order of 
magnitude of the error is about (^)^ — 10""^. 

A provision must be made for the initial conditions on the mode functions solution 
of (|3.2|, |3.3|) in the broken symmetry case. In this case the tree level squared mass is 



negative m| = — //^ and for initial values of the expectation value 0^(0) < 2/i^/A the initial 
configuration is below the classical "spinodal" and imaginary frequencies lead to the spinodal 
instabilities. We are not interested in this article in studying the time evolution of these 
instabilities but on the issues of non-linear relaxation. There are two ways to avoid the 
complex frequencies associated with these instabilities in the initial conditions: (i) one can 
choose an initial value 0^(0) > 2/i^/A or (ii) one can impose that the initial frequencies 
correspond to a positive mass term. This choice corresponds to preparing an initial gaussian 
density matrix of the modes with this given mass, under time evolution this packet spreads 
in function space and the time dependence of the width reflects the (quantum and thermal) 
fluctuations. In our numerical analysis we chose the later possibility with the positive mass 
squared given by the absolute value of the (negative) mass squared in the Lagrangian. 
We introduce now the following dimensionless variables: 



t]{t) 



A 



6 |m| 
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In terms of which the evolution equation ( |3.5| ) becomes 



q dq 



|f/.(r)P 



1 9 2 - lK(r)P - 1 



+ ^vir)l-q^dq-^^ 



Wn 



21 



A 



In 



A 



(3.6) 



where we chose the renormahzation scale k = \mq,\. 



1. Particle Production 



As in any time dependent situation, the concept of particle is ambiguous and has to be 
specified with respect to some particular state. We choose that state as the equilibrium 
ensemble at the initial time t = 0. The initial condition on the mode functions for the 
fluctuations 



^ |3.3| ) naturally determine the set of positive and negative energy states at 
this initial time. At this time the fluctuation operators may be expanded in this basis. The 
Fourier components of the fluctuation operators are thus written as 



1 



am - 4(0) 

%(0) - 4(0) 



The number operators at any time t are 

Tr4(0)a,^(0)p(t) _ Tr 4(t)a,Kt)p(0) 



Trp(O) 



Tr p(0) 
a^.it) = Uit)amU~\t) , 

Tr4(0)6,K0)p(t) _Tr4(t)6,,(t)p(0) 



Trp(O) 



Trp(O) 



bj:it) = U{t)bmU-\t) 



with U{t) the unitary time evolution operator. Following the arguments presented in refer- 
ence []T3[ we flnd that the creation and annihilation operators at time t are related to those 
at the initial time t = by a Bogoliubov transformation. In terms of the dimensionless 
variables introduced above we flnd 



Ar^,,(r) = (2^-l)Ar^,,(0) + (^-l) 



1 
2 



\V,ir)\' 
w',\V,{t)\^ 



+ 



where N^^rj{0) are the occupation numbers at the initial time t = and in the case given by 
the Bose-Einstein distribution functions, derivatives are with respect to the dimensionless 
variable r. Since in this section we will be working at zero temperature, only the last term 
will contribute to particle production. This term is recognized as the induced contribution. 
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It can be seen from the initial conditions on the wave functions that the induced contribution 
vanishes at the initial time. 

We will compute numerically the number of particles produced in a correlation volume 



X\t) = / d^kNk{t)/\m^\ 



(3.7) 



2. Energy Conservation: 

To one loop order the expectation value of the Hamiltonian (divided by the volume) 
gives the mean energy density 



E 



1 



W^{t) = k^ + M\t) 
wl{t) = k^ + m^{t) 



2 J (27r)3 [2WkiO) 
Vkf + wl{t)\Vkf^ 



\Uk\^ + W^it)\Uk 



(3.8) 

(3.9) 
(3.10) 



Using the one-loop equations for the zero mode and the mode functions it is straightfor- 
ward to find that the above energy density is conserved in time. Thus this (resummed) one- 
loop set of self-consistent equations is energy conserving. This is an important consistency 
check of the set of equations that guarantees that the numerical results are trustworthy. 

B. Inflaton Coupled to Fermions Only 

The fermionic Green's functions are constructed from the solutions to the homogeneous 
Dirac equation in the presence of the background field. The main ingredients and treatment 



is similar to that of fermions in presence of an electric field studied by Kluger et. al. p6 
Writing the independent solutions as 



W(i'2)(f,t) 



ifc-xrr(l,2) 
1^1. 



the mode functions obey 



^7o^ - 7- k-m^{t) 

d ^ 
«7o^ + 7 ■ k-m^[t) 



It is convenient to write the spinors as 
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«7o^ - 'J ■ k + m^{t) 



with u^^''^\ the spinor eigenstates of 70 with eigenvalues +1, —1 respectively. The 

functions fkit), gkif) obey the second order equations 



— + fc^ + m^(t) - im^it) 

— + /c^ + m^(t) + im^{t) 



fk{t) = , 
gk{t) = . 



We now need to append initial conditions. We will consider the situation in which the 
system was in equilibrium at time t < with the expectation value of scalar field being 0(0) 
and 0(0) = 0. Thus the fermion mass is constant and given by m^(0). We can now impose 
the condition that the modes fkit), Qkit) describe positive and negative frequency solutions 
for t < and, normalizing the spinor solutions to the Dirac equation to unity, we impose 
the following initial conditions 



fk{t < 0) 
Qkit < 0) 



-ikot 



2ko{ko + m40)) 
2/co(/co + m^(0)) 



ko = ^Jk^ + ml{0) . 

Equations ( |3.11| ) with these boundary conditions imply that 

9kit) = flit) . 



(3.11) 



The necessary ingredients for the zero temperature fermionic Green's functions are the fol- 
lowing 

s>{tx) = -^ E 
0=1,2 

a=l,2 

In the standard Dirac representation for the 7 matrices, we find 

•l+7o 



S>{t,t') = -if,{t)f;{t') [Wk{t)jo-7-k + m4t) 
S<it,t') = -tfl{t)fk{t') [Wfc*(t)7o-7-^-m^W 

■ hit) 



2 

1 - 7o 



>V,*(07o-7-fc + m^(0 
Wk{t')^o-l-k-m^{t') 



fkit) 
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With these ingredients the non-equihbrium fermionic Green's functions can be constructed 
as in section II. It is an important check that the equality TrS'^^(t,t) = Tt S^~{t,t) = 
TTSj^(t,t) = TrS^itjt) is fulfilled where the trace is over Dirac indices. It is also an 
important property that 

+ m4t)f;{t)) (tUt) + m^{t)h{t)) + efl{t)h{t) = 1 . (3.12) 

This is a consequence of the conservation of probability and may be checked explicitly. 
The evolution equation becomes 

0(t)+m|0(t) + ^0=^(t) + ^0(t) <x'(x,t) > -yTrS<(f,t;f,t) = 0, 
D 2 



Tr^<(f,t;f,t) = 2 



(2^ 



2k'\m\' 



(3.13) 



The expectation value of the Hamiltonian can again be computed in this (resummed) 
one-loop approximation. Using the equations of motion for the expectation value of the 
scalar field and those of the mode functions of the fermions, it can be shown to be conserved 
similarly to the bosonic case. 

The integral in ( ^.13 ) is divergent. The divergence structure may be understood in two 
different ways: carrying out a WKB expansion of the modes, just as was done in the bosonic 
case (this is a tedious and lengthy exercise), or alternatively calculating the closed loop in the 
presence of the background field. We carried out both methods using an upper momentum 
cutoff in the integrals and found 
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d^k 
(2^ 



i-2k'\m\' 



y 

"27r2 



^ 'm4t)+2ml{t)]\n- + 



47r2 



finite 



(3.14) 



where k is an arbitrary renormalization scale which will be chosen again to be k = m$ for 
numerical convenience. We will concentrate on the "chiral limit" in which the bare mass 
term for the fermions vanishes and m^{t) = y(l){t). In this limit the discrete symmetry 
(f) —>■ —(f) is explicit. The divergent terms in ( |3.14D are then identified as: mass renormaliza- 
tion (of (f)) (first term), wave function renormalization (second term) and coupling constant 
renormalization (third term). From now on, all quantities will be renormalized, and we drop 
the subscript R to avoid cluttering. For simplicity in the numerical analysis and in order to 
isolate the fiuctuations from the fermionic contribution from those of the scalar sector, we 
now neglect the contribution from the scalar loop. The scalar fiuctuations had already been 
analyzed in the previous section and in [|l^. The final renormalized equation of motion is 
now 



0(t) + ml (Pit) + ^(P^t) 



2y 



r 



47^2 



4>{t) + 2y' (P%t) 



d^k 
(2^ 
A 1 



In- 



. 



A^ 0(t) 
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With the purpose of a numerical analysis of this equations, it proves convenient to introduce 
the following dimensionless of variables 
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r]^y(f)/m^, r = nhi,f . q = k/m^ , g' ^ Xn/Qy'^ , g^y'^/ir'^, 
u,ir) = Mt)^u;t{u;t + y'nO)) , 
in terms of which the equations of motion are 
cfr] 



575 + 5' + ''' 



Jo 



A/m$ 



1 

2 



irj 

dq 

A V 



Ugir) = , Ug{0) = -j= , uq{0) = -i 



q^ + r)^{0) 



V2 



1 - 



1 \Ug[^) 



?7 + ^ (17 + 2 r/^) In 



A 



(3.15) 



Although the equations are finite when the ultraviolet cutoff is taken to infinity and can be 

written without the introduction of the cutoff by subtracting the integral with a lower limit 
cutoff given by the rcnormalization scale, numerically these integrals will have to be done 
by introducing an upper momentum cutoff anyways. Therefore we keep this UV cutoff in 
the equations. 



1. Particle Production 



The study of particle production is similar to the scalar case, with the only complication 
being the spinorial structure of the fermionic fields. At the initial time t = the quantized 
fermion operator is 



V'(f,0) 



d^k 
(2^ 



a=l,2 



Ak-x 



(3.16) 



in terms of the mode functions determined above. The number of fermions and antifermions 
are defined as 

[p(t)6<"»(0)6<°'(0)l 



Tr 



d^k 
(2^ 



E 



Tr 



TVp(O) 

p(t)4"^Ho)4")(o) 



r d^k ^ 



Trp(O) 



(2ny 

dH: 
{2ny 



E 



Tr 



TVp(O) 

p(o)4"^^w4"Ht) 



Trp(O) 



The time dependent coefficients are obtained by projecting the time dependent spinor so- 
lutions onto the positive and negative energy solutions at i = 0, and are related to the 
coefficients at t — via a Bogoliubov transformation 



E[^ 



(a) 

(/3)fcH 



(*)i.f 
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The identities 



/3 

ensure that the transformations preserve the anticommutation relations as they must. It 
is also a matter of algebra using the equations of motion for the mode functions to prove 
that the number of fermions minus antifermions is conserved for each k mode. In what 
follows, we only consider the number of fermions produced since fermions and antifermions 
are created in pairs. 

With some algebra, the Bogohubov coefficients can be expressed as the function of fkit) 
and fl{t). In terms of the dimensionless variables defined above we obtain the number of 
fermions within a correlation volume M^{t) = {T)/\m<j,\^ 



1 



27r2 
1 

4^ 



dqUHr) 



dq 



d 



q^ + ri^{0) -ri{0) 



d I 

+ ^(^) - V + ^^(0) - r7(0) 



(3.17) 



IV. NUMERICAL ANALYSIS 

The numerical analysis was carried out with a fourth order Runge-Kutta method for the 
differential equations and a 5-point Bode rule integrator for the fc-integrals. The typical step 
size in time was 1 — 2.10"^, and the typical step size in k was 10^^. The cutoff A/|m<j,| was 
varied between 75 and 200; we found no sensitivity to the cutoff in the range of parameters 
that we tested (see below). The code was very stable within the range of parameters tested. 



A. Inflaton Coupled to Scalars: Unbroken Symmetry Case 

Since the contribution of the one- loop quantum fluctuations of the inflaton $ have already 
been studied previously [jl3l and we want to study the contribution of the lighter fleld a, we 
will neglect the contribution from the self-interaction in the evolution equation. Thus we 
study the following equations obtained from equations (|3.3| , |3.6| ) in terms of dimensionless 
variables 

r)(r) + r^{r) + r^\r) + r^(r) q^dq 
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m<i, 



0, 



\ 



m 



6g 



A 



Notice that the factor (3/2) multiplying the logarithm in eq. ( p.6| ) and missing from ( [4.1| ) 
arises from the renormalization of the $-scalar loop that is not taken into account in ( [4.1D . 

Figures (l.a-c) show the unbroken symmetry case with the quantum fluctuations from 
the scalar loop of the a particles for the values y = 0; A/Svr^ = 0.2; g = X; = 
0.2 m$; r]{0) = 1.0; r]{0) = 0. Figure (l.a) shows rj^r) vs r, figure (l.b) shows A/'o-(r) vs r 
and figure (l.c) shows Afq^a- vs g for r = 120, similar graphs were obtained with snapshots 
at different (earlier) times. 

Figure (l.a) shows a very rapid, non-exponential damping within few oscillations of the 
expectation value and a saturation effect when the amplitude of the oscillation is rather 
small (about 0.1 in this case), the amplitude remains almost constant at the latest times 
tested. Figure (l.a) and figure (l.b) clearly show that the time scale for dissipation (from 
figure (l.a) is that for which the particle production mechanism is more efficient (figure 
(l.b)). Notice that the total number of particles produced rises on the same time scale as 
that of damping in figure (l.a) and eventually when the expectation value oscillates with 
(almost) constant amplitude the average number of particles produced remains constant. 
These figures clearly show that damping is a consequence of particle production. At times 
larger than about 40 mj^ (for the initial values and couplings chosen) there is no appreciable 
damping. The amplitude is rather small and particle production has practically shut off. If 
we had used the classical evolution of the expectation value in the mode equations, particle 
production would not shut off (parametric resonant amplification), and thus we clearly see 
the dramatic effects of the inclusion of the back reaction. 

In this unbroken symmetry case linear relaxation predicts a slow t^^ power law decay to 
an asymptotic finite amplitude because one particle decay is kinematically forbidden: the self 
energy contribution is a two-loop effect with one $ and two a particle cut and kinematically 
there is a one-particle pole below the three particle threshold. A slow power law linear 
relaxation asymptotically cannot be ruled out numerically because we have not continued 
the integration for longer times but clearly asymptotically the numerical result is compatible 
with linear relaxation. Figure (l.c) shows the distribution of particles created at the latest 
time T = 120 as a function of wave vector. Similar graphs were obtained with snapshots 
at different earlier times. The distribution is clearly non-thermal and skewed towards small 
momentum (in units of m$). These figures are qualitatively similar to those obtained in 
the self-interacting case in 



The asymptotic behavior is that of undamped oscillations 
of small amplitude. This is compatible with the result from the linear relaxation analysis 
because there is a one-particle pole below the three particle threshold resulting in undamped 
oscillations at large times. Linear relaxation predicts qualitatively the same results for the 
self-interacting scalar case and this case in the unbroken phase. The numerical results are 
consistent with this prediction. 
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However, for large amplitudes non-linear relaxation via particle production is very effec- 
tive and dramatically different from linear relaxation. 



In this case, we take m| = — |m||. 

As in the unbroken symmetry case, we only study the effect of the lighter a fluctuations. 
In the broken symmetry case, linear relaxation predicts an open decay channel for the 
inflaton, resulting in exponential relaxation for a long time and eventually relaxation with 
a power law. Therefore we should not expect a constant amplitude asymptotically biit an 
amplitude that eventually should relax to zero. In this case the renormalized equations for 
evolution are 



Figure (2.a-c) show ?7(r) vs r, A/'o-(t) vs r and Mq^aij = 200) vs q respectively, for A/Stt^ = 
0.2; g/X = 0.05; = 0.2 |m<i.|; 7^(0) = 0.6; r]{0) = 0. Notice that the mass for the 

linearized perturbations of the $ field at the broken symmetry ground state is \/2 |m$| > 
2mo-. Therefore, for the values used in the numerical analysis, the two-particle decay channel 
is open for linear relaxation. For these values of the parameters, linear relaxation predicts 
exponential decay with a time scale Trei ~ 300 (in the units used). Figure (2. a) shows very 
rapid non-exponential damping on time scales about six times shorter than that predicted 
by linear relaxation. The expectation value reaches very rapidly a small amplitude regime, 
once this happens its amplitude relaxes very slowly. Within our computing time limitations 
we could not confirm that there is exponential relaxation in the small amplitude regime (for 
r > 100) but clearly there is a striking difference with the unbroken symmetry case. The 
influence of open channels is evident, however in the non-linear regime relaxation is clearly 
not exponential but extremely fast. Although we cannot confirm the exponential (or power 
law) relaxation numerically in the small amplitude regime, the amplitude at long times seems 
to relax to the expected value, shifted slightly from the minimum of the tree level potential 
at ?7 = 1. This is as expected from the fact that there are quantum fluctuations. Figure 
(2.b) shows that particle production occurs during the time scale for which dissipation is 
most effective, giving direct proof that dissipation is a consequence of particle production. 
Asymptotically, when the amplitude of the expectation value is small, particle production 
shuts off. We point out again that this is a consequence of the back-reaction in the evolution 
equations. Without this back-reaction, as argued above, particle production would continue 
without indefinitely. Figure (2.c) shows that the distribution of produced particles is very 



B. Inflaton Coupled to Scalars: Broken Symmetry Case 
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far from thermal and concentrated at low momentum modes k < |m$|. This distribution is 
qualitatively similar to that in the unbroken symmetry case, and points out that the excited 
state obtained asymptotically is far from thermal. 



C. Inflaton Coupled to Fermions Only: Unbroken Symmetry Case 

Here we treat the case y ^ 0; g = 0. 

The renormalized evolution and particle production equations in this case are given by 
eq. ( p.l5| ) and eq. (|3.17| ), respectively, in terms of dimensionless quantities. Figures (3.a-c) 
show ?7(r) vs. r, {t) vs r and A/J(r = 200) respectively for the values of the parameters 
~ 0; Z/^/tt^ = 0.5; = 1.0; //(O) = 0.6; ?7(0) = 0. One observes from these 

figures that after a rather brief period of initial damping of just a few oscillations of the scalar 
field, the dissipative mechanism shuts-off. Figure (3.b) shows that during this time scale 
fermion-antifermion pairs are being produced but then the number of produced particles 
saturates and oscillates with a small amplitude. Figure (3.c) shows that the distribution of 
fermions produced is peaked at very low momentum {k < mq>) and with a maximum value 
of 2, which is the total number of degrees of freedom per A;- wave vector. These numerical 
results expose the physics of Pauli blocking very clearly; the available low momentum modes 
are occupied and no more fermion-antifermion pairs can be produced. Pauli blocking shuts 
off particle production and dissipation very early on. We have obtained snapshots of the 
particle number as a function of momentum for different times, and they all present the 
same picture. 

This result is markedly different from the prediction of linear relaxation. At zero temper- 
ature, eq. ( |2.48| ) for linear relaxation predicts exponential damping with a (dimensionless) 



time scale Trei ~ 10.6 for the values of the parameters chosen above. The difference between 
the non-linear evolution and that predicted by linear relaxation is explained by the Pauli 
blocking phenomenon. Very early in the evolution, the low momentum available fermionic 
states were filled with produced fermions. Once these states have filled, damping and par- 
ticle production shuts off. This Pauli blocking effect is explicit in eq. (|2.48| ) but there it 
appears from finite temperature effects. In the non-linear relaxation case, we began at zero 
temperature, but an excited state quickly ensues because of particle production. 

This analysis reveals that for large amplitudes of the scalar field, fermions will be rather 
ineffective in dissipation and damping because of Pauli blocking even at zero temperature. 
The time scales for dissipation obtained from the fermion self-energies, which apply to the 
case of linear relaxation are completely unrelated to the time scales for non-linear dissipative 
processes, even asymptotically, because the fermionic states are Pauli blocked. 



D. Thermalization: the argument 



The self-consistent resummed one-loop approximation used in our analysis does not in- 
corporate collisional processes. These will appear in the next (two loop order). The distri- 
bution of the particles produced is clearly very far from thermal, with typically very many 
particles in the low momentum modes k < m^. Although a full study of thermalization will 
necessarily involve the two loop contribution (and consequently memory effects that will 
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ultimately result in an extensive numerical effort), we can argue on the basis of an effective 
Boltzmann description. We can take the distril:)ution of particles obtained from the numer- 
ical analysis at long times (times longer than the non-linear relaxation scale) as an input in 
a Boltzmann equation. The integration of this Boltzmann equation will determine the time 
scales for coUisional thermalization. Because the couplings will be typically very small and 
the initial occupations are rather far from thermal equilibrium these collisional time scales 
will be very long (as compared to the non-linear relaxation scales). Collisional relaxation 
does not change the number of particles and is energy conserving (in flat Minkowsky space), 
just re-distributes the occupation. Therefore a thermal equilibrium state will result with 
the total number of particles given by the number of produced particles obtained during the 
fast period of induced amplification and the final reheating temperature will be determined 
by this number of particles. Including the collisional relaxation via the Boltzmann equation 
during the period of induced amplification will not modify the time scales substantially. The 
reason is that initially the occupation numbers are rather small and Boltzmann processes are 
suppressed. Since the rate of particle production depends on the couplings and the initial 
amplitudes, for large values of the amplitudes the rate of particle production will be much 
larger than the rate of equilibration via Boltzmann processes (second order in the couplings 
and independent of the initial amplitude). Thus collisional (Boltzmann) processes will be- 
come relevant for thermalization at times much longer (for weak coupling) than the times 
during which most of the particle production takes place. If the rate of thermalization ob- 
tained from this Boltzmann equation is much larger than the expansion rate of the universe, 
then the reheating temperature can just be estimated from the total number of particles 
produced during the period of induced amplification. However if this collisional thermaliza- 
tion rate is much smaller than the expansion rate there will be substantial redshifting of the 
temperature and a full numerical analysis will have to be carried out to obtain a reliable 
estimate of the reheating temperature. 



With the arguments on thermalization given above in mind, and assuming that the ther- 
malization rate is much larger than the expansion rate (no red-shifting of the temperature) , 
we can use the numerical analysis presented above to obtain an estimate of the reheating 
temperature. 

The important conclusion of the numerical integration is that the density of particles 
produced during the period of induced amplification is 



where is a mass in between the infiaton and the a field masses, c a constant of (9(1 — 10). 
This constant depends on the product of the coupling constants and the initial amplitudes as 
can be seen from the mode equations. It also depends approximately logarithmically on the 
coupling constants because of the time scales for the shutting-off of the induced amplification 
mechanism. This can be seen as follows: from the expression for the number of (bosonic) 
particles produced we find the rate of particle production (for example for the a field) 



E. Estimate of the Reheating Temperature 
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This expression clearly shows that the rate of particle production depends on the produc 
of the coupling and the initial amplitude (squared for the interaction considered), since the 
mode functions also depend on this combination. The total number of particles is given by 
this rate integrated up to the time at which the amplification mechanism is shut-off by the 
back reaction. Un upper bound to this time is estimated as follows: in the early stages, 
the mode functions increase almost exponentially (because of the resonances) this growth 
is modified by the damping of the amplitude of the zero mode. Since we find that the back 
reaction shuts-off the particle production when the one loop term becomes of the order of the 
tree-level term, this gives a time scale that is logarithmic in the coupling. Thus even for very 
small couplings but very large initial amplitude of the inflaton zero mode (and this is the 
situation in cosmological scenarios) the constant c will be of 0{1 — 10) and can be enhanced 
by giving the inflaton a larger initial amplitude. Under the assumption of thermalization on 
time scales shorter than the expansion scale and that the mass of scalar s and fermions are 
much smaller than the inflaton mass an estimate of the reheating temperature is obtained 
as 

Treh ^ AT^ ^ TUa. (4.3) 



V. CONCLUSIONS AND IMPLICATIONS 

The reheating and thermalization processes in inflationary universe models occur very 
far from equilibrium and must be studied in their full complexity, eventually numerically. 
The methods developed within real-time non-equilibrium field theory allow a consistent 
treatment that we used to obtain the equations of evolution for the expectation value of the 
inflaton field. These equations are non-perturbative and take into account the back reaction 
effect of both the self couplings of the inflaton, as well as the couplings to lighter scalars and 
fermions. 

We have examined these evolution equations, both in the linear regime, in which the 
amplitude of the field is small and in the opposite extreme, the non-linear regime, which is 
intrinsically non-perturbative. 

In the linear regime, we have seen that if the inflaton mass is above the two particle 
threshold, there is some damping behavior for a period of time, and indeed, this damping 
is governed by the total width of the inflaton. However, at late times the behavior is 
dominated by power law decay rather than exponential damping. Furthermore, we have 
noted the important distinction between the time scale for the relaxation of the expectation 
value of the inflaton as opposed to thermalization of the produced particles. It is clear that 
as treated in this work, relaxation can occur well before thermalization, i.e. interactions that 
drive the particle distributions towards a thermal one, has a chance to become relevant. 

The more interesting case, and the one which would most likely be relevant to the 
reheating problem in models such as chaotic inflation 0, is that of when the inflaton is in 
the non-linear regime. Unlike the elementary approach to reheating, in which the "decay" of 
the coherent inflaton oscillations occured a time tdecay ~ T"^, where F is the inflaton decay 
width, after the end of the slow-roll regime, particle production in the non-linear regime 
begins at very early times and then shuts off. Indeed, by the time the inflaton expectation 
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value reaches its asymptotic state of oscillations around its minimum, the period of particle 
production is essentially over. This will be the stage of thermalization via collisions, however 
we find that the spectrum of produced particles is extremely non-thermal, and collisional 
relaxation towards an equilibrium thermal state may take a long time. 

An inflationary model that cannot reheat the universe to at least nucleosynthesis tem- 
peratures is wrong. Thus it is clearly important to accurately compute the reheating tem- 
perature in inflationary models. What our analysis here shows is that this computation is 
of necessity much more complicated than has previously been thought but perhaps more 
importantly that there are different time scales. Although so far we have studied the situa- 
tion in flat Minkowski space, recent investigations by Kaiser reveal that these non-linear 
processes will still be very effective in inflationary cosmologies. 

The first stage which occurs rather fast is that of particle production and damping of 
oscillations as a result of induced amplification. The second stage, that begins typically 
when the amplitude of the inflaton is rather small and it oscillates with almost constant 
amplitude at the bottom of the potential is that of one-particle decay and of thermalization 
via collisional relaxation. 

The thermalization time will have to be studied setting up a (quantum) Boltzmann 
equation using the distributions of produced particles at the end of the induced amplification 
stage as input for the Boltzmann evolution. Since these initial distributions are very far 
from equilibrium, it may take many collisions to relax to a thermal equilibrium state. If 
the couplings are very small (and this will clearly depend on the models) the thermalization 
time may be very large. 

Following the particle numbers to the point at which they become thermal would then 
allow us to pick off the final reheating temperature. In the weak coupling regime, as men- 
tioned above, the thermalization rate may be smaller than the expansion rate, allowing 
for significant redshifting of the total energy density, leading to a low reheat temperature, 
perhaps of the order of the inflaton mass or less. 

The particular case of large thermalization rates compared to the expansion rate, al- 
lows for a quick estimate of the reheating temperature. In this case we can assume that 
thermalization occurs without red-shifting of the energy which is then conserved during the 
thermalization time. By taking the energy density at the beginning of this stage to be aT^ 
with a the Stephan Boltzmann accounting for the particle statistics one can then obtain an 
estimate of the reheating temperature, but the justification for this should ultimately arise 
from a deeper understanding of the time scales involved. 

Our main point in all of this, though, is that self-consistent calculations taking into 
account the non-linearity of the dynamics are necessary in order to extract the reheating 
temperature in a reliable way. Our work clarifies what the relevant particle production 
mechanisms are in these theories and provide a consistent and implementable framework of 
calculation. We are now extending these studies to FRW cosmologies. 

We emphasized in this study that thermalization and particle production are fundamen- 
tally different processes and in the non-linear regime likely to occur on widely different time 
scales. 
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A. A possible solution to the Polonyi problem: 



As mentioned before a problem of much current interest in which this work may have 
some imphcations is the "post-modern Polonyi problem" concerning flat directions for some 
of the moduli fields in string theories In these scenarios the Polonyi fields are very 

light (of order of a Tev) with couplings to the particle physics sector that is extremely small 
(order of 10~^^). A calculation of the reheating temperature based on the "old" scenario of 
one particle decay give T^eh ~ 0.01 Mev , well below the nucleosynthesis scale. 

This may not necessarily be true given our analysis. Since particle production occurs 
not just during oscillations around the minimum now but throughout the evolution of the 
the moduli field. Even for couplings 0{m^/Mpi) since the initial amplitudes of these moduli 
field are 0{Mpi) then the mechanism of particle production via induced amplification will 
be very efficient and our analysis leads to an estimate of the reheating temperature Treh ~ 

~ 1 Tev that is larger than the electroweak scale, alleviating the Polonyi problem. Thus 
this mechanism may very well provide a solution to the Polonyi problem reconciling moduli 
field cosmology with Standard Big Bang cosmology. 
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APPENDIX A: A PEDAGOGICAL EXERCISE 

In this appendix we show explicitly how to implement the tadpole method for the case of 
the infiaton coupled to a lighter scalar field via a simple trilinear coupling. We will not worry 
about renormalization issues here, though they were, of course, dealt with in the text. The 
formalism for non-equilibrium quantum field theory has already been described several times 
in the literature [^. A path integral representation involves an integration along a path in 
the complex time plane with forward, backward branches and if the initial density matrix 
was that of an equilibrium system at an (initial) temperature T also an imaginary time 
branch. For real time correlation functions the imaginary time branch does not contribute 
and only determines the boundary conditions on the Green's functions. 

We use the tadpole method to study the time evolution of the expectation value 

of the infiaton. 

' ' Trp(O) Trp(O) ' ^ ^ 
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where are the fields defined on the forward and backward branches respectively. We set 



(A2) 



where ^ire field fiuctuation operators defined along the respective branches. 

The effective evolution equation for the background field t) follows from the condi- 
tion. 



< x^{x,t) >= . 



(A3) 



Treating the linear and non-linear terms in as interactions and imposing the condition 
( |A3| ) consistently in a perturbative or loop expansion, one obtains expressions of the form 
(here we quote the equation obtained from < x+(x, t) >= ) 



dx'dt' 



< X+(x,t)x+(x ,0 > 0++(x',t')+ < X^{x,t)x-ix,t') > 0+-{x',t')] = (A4) 



and similarly for < x^(af, t) >= 0. The are in general integro-differential operators 
acting on the background field. 

Because the Green's functions < x"'"(x, t)x^{x', t') >, etc. are all independent one obtains 
the equations of motion in the form 



C)++(f' = , 0+-{x',t') =0 , 0^+{x',t') = , O—{x',t') = 0. 



(A5) 



It is a consequence of the properties of the non-equilibrium Green's function (see below) 
and ultimately a consequence of unitarity that all the integro-differential operators O are 
the same. 

Consider the Lagrangian density 



Ci<^,a) = Com + Coia) + git) a' 



(A6) 



with the Cq being the free field Lagrangian density (with respective mass terms) and 
have allowed the coupling to depend on time. The non-equilibrium path integral requires 
£($"'", o""*") — £($~,o"~). In the tadpole method we write $^(a;, t) = x^(2^, t) + (pit) and 
identify 0(t) as the (non-equilibrium) expectation value of the field $ . This identification 
then requires that < xi^^t) >= where the expectation value is in the non-equilibrium 
density matrix with the path integral representation along the contour in complex time. 
After this shift, the action reads: 



L = J £xdt {Coix^) + Coia^ 

[x'^^fy^ ^ )} • 



-0 - m|0l + (7(t)0(t)(a+)2 + (7(t)x+(fx^ 



The linear term in is included as a perturbation. The first contribution to the equation 
of motion is obtained from this linear term, from the condition < x^i^i t) >= obtains 
to this order 



dt' {< X+(t)x+(t') > i^) [-4> - <<!>] - < X^it)x-it') > i^) [~4> " ml<p\ } = , (A7) 
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where we have suppressed the spatial arguments. Because the correlation functions 

< X'^{'t)x^{'t') >) < X^{'t)x^it') > si-re independent, one obtains the tree level equa- 
tions of motion. It is straightforward to see that the same is obtained by imposing 

< x~(a;, t) >= 0. In an amplitude expansion (an expansion in powers of (j){t)) the one 
loop correction to the equation of motion is obtained by expanding (the exponential of) 



to first and second order. The first order gives a tadpole 



contribution, the second order needs one vertex with x, the other with (p. One obtains 



dt' < xHt)x^{t') > {(^) (-0 - + m')) < {^Ht')f 



> + 



< a 



-{ma+{t")f 



> - < 



J dt"{ig{t")f 

J dt' < xHt)x~{t') > {{^) (-0 - mU) + {ig{t')) < {a-f{t') > - 

J dt"{-ig{t")f [< {a-{t')na-{t")f > 0(O- < {a-{t')na^{t")f >] m] = . 

The expectation values are computed using Wick's theorem and using the free-field Green's 
functions of section II. The tadpole (time independent) is absorbed in a shift of the expec- 
tation value. The coefficient of < x^x^ >, < X^X > must vanish independently because 
these Green's functions are independent and must vanish at all times. From the Green's 
functions (and more generally from the formal time contour integral) it is seen that the 
equations obtained are identical. If the expectation value is translational invariance, the 
spatial integrals set the momentum transfer to zero. For example using the zero tempera- 
ture Green's functions of section I, one finds that the term that has the non-local (in time) 
correlation functions (last term) becomes 

dt"{tg{t")f [< {a+{t')f{a^{t")f >-< {a+{t')na-{t")f >] <P{t") = 

d'^k 



: / dt"{ig{t")f<t>{t")e{t' 



t" 



For the resummed one-loop approximation the term gd>{t){a'^Y is absorbed in the mass 
term for the a field and only the xc^ terms are considered as perturbation. In this case 
the 1-loop contribution is obtained from the tadpole term < (a'^{t)Y >, which is now time 
dependent and obtained from the non-equilibrium Green's functions constructed with the 
mode functions for the time-dependent mass as in section IV. One now finds the evolution 
equations 



(f){t) + ml (f){t) + g{t) < a^{x, t)>^0, 



d^k 
Vk{t) 



<a\x,t) >= J ^sGt;At,t) = J 

0, 



(27r)3 2wk 



coth 



Vk{0) = 1 , 14(0) = -iwk = -i^k^ + m-'iO) , 
m^(t) ^ma + g(t)(l)(t) . 
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Figure Captions: 

Fig.(l.a) Scalars: unbroken symmetry ?](r) vs r for the values of the parameters 
y = 0; A/87r2 = 0.2; ^/A = 1; m„ = 0.2 m^; 77(0) = 1.0: ?7(0) = 0. 

Fig.(l.b): Nair) vs. r for the same value of the parameters as figure (l.a). 

Fig.(l.c): Mq,a{T — 120) vs. q for the same values as in figure (l.a). 

Fig. (2. a) Scaleirs: broken symmetry ?7(r) vs r for the values of the parameters y — 
0; A/Stt^ = 0.2; g/\ = 0.05; m„ = 0.2 \m^\; r/(0) = 0.6; 77(0) = 0. 

Fig.(2.b): A/"o-(r) vs. r for the same value of the parameters as figure (2. a). 

Fig.(2.c): Mq^air = 200) vs. q for the same values as in figure (2. a). 

Fig. (3. a) Fermions: unbroken symmetry 7]{t) vs t for the values of the parameters 
^ = 0; yVTT^ = 0.5; X/6y^ = 1; = 0; 7^(0) = 1.0; 77(0) = 0. 

Fig.(3.b): {t) vs. r for the same value of the parameters as figure (3. a). 

Fig.(3.c): N^{t — 200) vs. q for the same values as in figure (3. a). 
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